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Abstract 



The shape of the transverse momentum distribution of W bosons (pj?) 
produced in pp collisions at yfs = 1.8 TeV is measured with the D0 detector 
at Fermilab. The result is compared to QCD perturbative and resumma- 
tion calculations over the pYf range from — 200 GeV/c. The shape of the 
distribution is consistent with the theoretical prediction. 
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The transverse momentum (py ) of W intermediate vector bosons produced in proton- 
antiproton collisions is due to the production of one or more gluons or quarks along with 
the boson. At low transverse momentum (p^ < 10 GeV/c), multiple soft gluon emission 
is expected to dominate the cross section. A soft gluon resummation technique [1-5] is 
therefore used to make QCD predictions. At high transverse momentum (pjf > 20 GeV/ c), 
the cross section is dominated by the radiation of a single parton with large transverse 
momentum. Perturbative QCD M is therefore expected to be reliable in this regime. A 
prescription |3J has been proposed for matching the low and high p^ regions to provide 
a continuous prediction for all p 1 ^ . Thus, a measurement of the transverse momentum 
distribution may be used to check the soft gluon resummation calculations in the low p^ 
range, and to test the perturbative QCD calculations at high p^ . 

The transverse momentum spectrum of W bosons has been measured previously by the 
UAl 0, UA2 ||, and CDF || collaborations, but with smaller data samples than the 
one used here. This paper presents a measurement of the shape of the pt spectrum of W 
bosons produced in pp collisions at y/s = 1.8 TeV with the D0 detector [|I(J at Fermilab, 
and extends the p^ range of the previous measurements. The data come from a sample of 
12.4 ± 0.7 pb^ 1 collected during the 1992-1993 run. A measurement of the inclusive cross 
section for W and Z boson production based on the same data set has been reported |TT[ 
and agrees with QCD predictions. 

This measurement uses the decay mode W — > ev. Electrons were detected in a hermetic 



uranium/liquid- argon calorimeter with an energy resolution of about 15%/JE(GeV). The 
calorimeter has a transverse granularity of At] x A0 = 0.1x0.1, where rj is the pseudorapidity 
and (f) is the azimuthal angle. Electrons were accepted in the central pseudorapidity region 
only, \rj\ < 1.1, to keep the background contamination from multijet events at a reasonably 
low level for high values of p^. . The transverse momentum of the neutrino was calculated 
using the calorimetric measurement of the missing transverse energy (Et) in the event. We 
take the pr to be the sum of the electron and neutrino transverse momenta, measuring it 
only from the recoiling hadrons. The analysis used a single electron trigger, which required 
one electron with transverse energy (Et) greater than 20 GeV. 

The offline electron identification requirements consisted of four criteria: (i) the electron 
had to deposit at least 95% of its energy in the electromagnetic calorimeter (21 radiation 
lengths deep); (ii) the transverse and longitudinal shower shapes had to be consistent with 
those expected for an electron ||12|| ; (iii) a good match had to exist between a reconstructed 
track in the drift chamber system and the shower position in the calorimeter; and (iv) the 
electron had to be isolated from other energy deposits in the calorimeter, with I < 0.1. This 
isolation variable is defined as / = [-&tot(0.4) — £ , EM (0.2)]/£ , EM ;(0.2), where E T q T (0A) is 
the total calorimeter energy inside a cone of radius y/ Arj 2 + A0 2 = 0.4 and £'em(0.2) is the 
electromagnetic energy inside a cone of radius 0.2. To select the W boson candidate sample, 
we required one electron with Et > 25 GeV, and E T > 25 GeV. Events having a second 
electron with E T > 20 GeV that satisfies criteria (i), (ii), and (iv) were excluded from the 
candidate sample as possible Z — > e + e~ events. Criterion (iii) was not applied to this second 
electron in order to allow for possible tracking inefficiencies. These selection criteria yielded 
7132 W — » ev candidates. 

The trigger and selection efficiencies were determined using Z — > e + e~ events in which 
one of the electrons satisfied the trigger and selection criteria. The second electron then 
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provided an unbiased sample with which to measure the efficiencies. No dependence of the 
trigger or selection efficiency on p^ was found, to an accuracy of 5%. 

A Monte Carlo program was used to simulate the D0 detector response and calculate 
the kinematic and geometric acceptance as a function of p™ ' . The detector resolutions used in 
the Monte Carlo were determined from data, and were parametrized as a function of energy 
and angle. The relative response of the hadronic and EM calorimeters was studied using the 
transverse momentum of the Z boson as measured by the pt of the two electrons compared 
to the hadronic recoil system in the Z event. This parametrized representation of the 
D0 detector was used to smear the theoretical prediction by detector effects and compare 
it to our measured p? . We prefer this method of comparison to a standard unfolding 



procedure |L4| that proved to be sensitive to the choice of the prior distribution function. 
This sensitivity is caused by the Jacobian zero in dN/dpf at the origin, which induces a 
peak that appears near p^ w 4 GeV/c after it is broadened by our p^ resolution. Only 
below 4 GeV/ c do the true p?f distributions predicted by the two available models |||| 
exhibit a difference, which is masked in the data by these same resolution effects. 

The dominant source of background in the W boson sample was multijet events where one 
or more of the jets fluctuated to fake an electron. Some multijet events also have significant 
$t due to fluctuations and mismeasurements of the jet energies. This could fake a neutrino 
from W boson decay. The amount of multijet background in the W — > ev candidate sample 
was determined by first defining a "loose" event sample which had the same selection criteria 
as the candidate sample except that electron identification criteria (i) and (ii) were not 
applied. This loose sample consisted of N s signal events and multijet background events. 
The W — > ev candidate sample (described above) consisted of e s N s signal events and e\,N h 
multijet background events, where e s and are the electron selection efficiencies for the 
candidate relative to the loose samples, for signal and background respectively. We obtained 
e s from the Z — ► e + e~ sample, and e\ } from events with < 15 GeV. The total multijet 
background was determined to be (4.2 ± 2.3)%. The shape of the multijet background as a 
function of p^ was determined by repeating the procedure in different p^ bins. 

To cross-check the multijet background estimate, the transverse mass and the $t dis- 
tributions of the final data sample were compared to a model of the expected distributions 
obtained from a combination of W — *> ev Monte Carlo events plus the estimated multijet 
background. The comparison was performed in three p^ bins: — 30 GeV/ c, 30 — 60 GeV/ c, 
and > 60 GeV/ c; the number ofW^ev candidates in each bin was 6726, 282, and 124. The 
amount of multijet background in each p^ bin was estimated as (2.9 ±1.6)%, (20.9 ±11.7)%, 
and (38.3 ±21.5)%, respectively. Figure [l] shows the results of the comparisons. The distinc- 
tive shape of the transverse mass distribution for multijet background arises from applying 
the kinematic cuts and the minimum p)^ requirement to a sample predominantly composed 
of dijet events that lie back-to-back in the transverse plane. The goodness of the fit be- 
tween the data and the model is evaluated by performing a x\ t es t [0- The numerical 
results for Xx/d.o.£. for each fit in Fig. | are (a) 20.2/25, (b) 26.4/19, (c) 4.7/7, (d) 5.4/9, 
(e) 3.0/13, and (f) 2.1/7. They correspond to fit probabilities of (a) 74%, (b) 12%, (c) 69%, 
(d) 80%, (e) 99%, and (f) 95% respectively. We therefore conclude that the tests show good 
agreement between the data and the expectation. 

The normalized multijet background was subtracted bin by bin from the W boson candi- 
date sample transverse momentum spectrum. Additional corrections (all less than 5%) were 
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made to account for top quark background events and for Z — > e + e~ events where one of the 
electrons was lost or not identified. Since p^ was measured from the recoiling hadrons, the 
events originating from W — > tv (where r — > ez/z/) contributed properly to the differential 
distribution; this source of background therefore was included in the Monte Carlo simulation 
of the p^ distribution. 

The normalized distribution of the W boson transverse momentum (j^-^w) is shown in 

Fig. |2] and given in Table |. The largest contributions to the systematic error in the p^ mea- 
surement are the uncertainty in the magnitude of the multijet background, the uncertainty 
in the hadronic recoil energy scale factor and resolution used in the detector simulation, and 
the uncertainty in the selection efficiency These are all added in quadrature since they are 
independent. We compare our experimental result to the theoretical prediction [3| computed 
using the MRSA' [HJ parton distribution function and smeared for detector resolutions. The 
measurement and the prediction are independently area-normalized to unity. These points 
are used to perform a more detailed comparison between data and theory by plotting the 
ratio (Data-Theory) /Theory, which is shown in Fig. |3|. Data and theory differ by less than 
a half of a standard deviation above p^ of 60 GeV/ c. We therefore conclude that the shape 
of the distribution is consistent with the theoretical prediction. 

In summary, we have measured the shape of the transverse momentum distribution of 
W bosons produced in pp collisions at y/s = 1.8 TeV, and have found that it is consistent 
with the combined QCD perturbative and resummation calculations. 

We thank the staffs at Fermilab and collaborating institutions for their contributions 
to this work, and acknowledge support from the Department of Energy and National Sci- 
ence Foundation (U.S.A.), Commissariat a L'Energie Atomique (France), State Commit- 
tee for Science and Technology and Ministry for Atomic Energy (Russia), CAPES and 
CNPq (Brazil), Departments of Atomic Energy and Science and Education (India), Col- 
ciencias (Colombia), CONACyT (Mexico), Ministry of Education and KOSEF (Korea), and 
CONICET and UBACyT (Argentina). 



6 



REFERENCES 



* Visitor from Universidad San Francisco de Quito, Quito, Ecuador. 
t Visitor from IHEP, Beijing, China. 

[1] J. Collins, D. Soper, and G. Sterman, Nucl. Phys. B250, 199 (1985); J. Collins and 
D. Soper, Nucl. Phys. B193, 381 (1981), Nucl. Phys. B197, 446 (1982), and erratum 
Nucl. Phys. B213, 545 (1983). 

[2] C. Davies, B. Webber, and W. J. Stirling, Nucl. Phys. B256, 413 (1985); C. Davies and 
W. J. Stirling, Nucl. Phys. B244, 337 (1984). 

[3] G. Altarelli, R. K. Ellis, M. Greco, and G. Martinelli, Nucl. Phys. B246, 12 (1984). 

[4] P. B. Arnold and R. Kauffman, Nucl. Phys. B349, 381 (1991). 

[5] G. A. Ladinsky and C. P. Yuan, Phys. Rev. D 50, 4239 (1994). 

[6] P. B. Arnold and M. H. Reno, Nucl. Phys. B319, 37 (1989); R. J. Gonsalves, J. 

Pawlowski, and C-F. Wai, Phys. Rev. D 40, 2245 (1989). 
[7] UA1 Collaboration, C. Albajar et al, Z. Phys. C 44, 15 (1989). 
[8] UA2 Collaboration, J. Alitti et al, Z. Phys. C 47, 523 (1990). 

[9] CDF Collaboration, F. Abe et al, Phys. Rev. Lett. 66, 2951 (1991), and Phys. Rev. 
Lett. 67, 2937 (1991). 
[10] D0 Collaboration, S. Abachi et al, Nucl. Instrum. Methods A 338, 185 (1994). 
[11] D0 Collaboration, S. Abachi et al, Phys. Rev. Lett. 75, 1456 (1995). 
[12] D0 Collaboration, S. Abachi et al, Phys. Rev. D 552, 4877 (1995). 
[13] D0 Collaboration, B. Abbott et al, Fermilab-Pub-97/328-E, |hep-ex/9710007| (October 



1997), submitted to Phys. Rev. D and Fermilab-Pub-97/422-E, |hep-ex/9712"029| (De- 
cember 1997), submitted to Phys. Rev. D. 
[14] G. D'Agostini, Nucl. Instrum. Methods A 362, 487 (1995). 
[15] S. Baker and R. D. Cousins, Nucl. Instrum. Methods A 221, 437 (1984). 
[16] A. D. Martin, R. G. Roberts, and W. J. Stirling, Phys. Rev. D 50, 6734 (1994); 51, 
4756 (1995). 



7 



TABLE I. The W boson transverse momentum distribution corresponding to Fig. 2. The 
column labeled "Stat Error" shows the statistical uncertainty; "Syst Error" shows the systematic 
uncertainty in background and efficiency; "Detector Error" shows the systematic uncertainty in 
the detector modelling; "Total Error" is the sum in quadrature of the previous three columns. 
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FIG. 1. The transverse mass (left) and $t (right) distributions for three pjf bins. The points 
are the D0 data. The solid histogram is the sum of the Monte Carlo signal and the estimated 
background. The dotted histogram is the estimated background alone. 
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FIG. 2. The W boson transverse momentum spectrum, showing the D0 result (solid points) 
with statistical uncertainty. The theoretical calculation by Arnold and Kauffman [4], smeared for 
detector resolutions, is shown as two lines corresponding to the ±la variations of the uncertainties 
in the detector modelling. Within each bin, the values are plotted at the mean pjf. The fractional 
systematic uncertainty on the data is shown as a band in the lower portion of the plot. The values 
of the uncertainties for different pj? bins are 100% correlated with each other. Upward fluctuations 
in the magnitude of the multijet background cause the widening observed in the band at about 60 
and 100 GeV/c in p$ . 
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FIG. 3. The ratio (Data-Theory)/Theory shown as a function of pjf with its statistical un- 
certainty as error bars. Within each bin, the values are plotted at the mean . The theory 
corresponds to reference [4], smeared for detector resolutions. The systematic uncertainties from 
data (background and efficiency) and from the detector modelling are added in quadrature and 
shown as a band. 
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